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Abstract One of the possible sources of gamma-ray bursts (grbs) are merging, compact neutronstar binaries. 
More than 90% of the binding energy of such a binary is released in the form of gravitational waves (gws) in the 
last few seconds of the spiral-in phase before the formation of a black hole. In this article we investigate whether a 
fraction of this GW energy is transferred to magnetohydrodynamic waves in the magnetized plasma wind around 
the binary. Using the 3+1 orthonormal tedrad formalism, we study the propagation of a monochromatic, plane 
fronted, linearly polarized GW perpendicular to the ambient magnetic field in an ultra-relativistic wind, first in 
the comoving and then in the observer frame. A closed set of general relativistic magnetohydrodynamic (grm) 
equations is derived in the form of conservation laws for electric charge, matter energy, momentum and magnetic 
energy densities. We linearize the GRM equations under the action of a monochromatic GW, which acts as a driver 
and find that fast magneto-acoustic waves grow, with amplitudes proportional to the GW amplitude and frequency 
and the strength of the background magnetic field. 
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1. Introduction 

The coupling between gravitational waves and electro- 
magnetic waves (emws) in a magnetized vacuum has 
been investigated ext ensively over the past 40 years by 
a number of authors ( Gertsenshteinlll96lt lLupa novlll967r, 
iBoccaletti et"aH Il97ft IZerdovichl Il973t Iberlachl 11974^ 
These studies demonstrate the coherent excitation of 
emws by a monochromatic GW propagating perpendic- 
ularly to a background magnetic field. 

The first calcula tions including the pr esence of a 
plasma were done bv lMacedo fe Nelson! I|1983|) . who found 
a coupli ng of GWs to ordinary and extraordinary emws, 
whereas |Brodin fe Marklundl l)l999(l derived the paramet- 
ric excitation of Langmuir wave s by a GW propagatin g 
through unm a gnetize d plasma. In lMarklund et all l|2000|l ; 
IBrodin et all fcOOlal). t he authors adopt the 3+1 te drad 
formalism l )Thorne et allll986tlEilis fc van Elslll998|) and 
show that the dispersion relation in a tenuous plasma dif- 
fers only from the vacuum solution by a small wavenum- 
ber shift. In a plasma, however, subsequent non-linear 
conversions such as harmonic generation, might allow 
the EMW energy to escape as radiation with frequencies 
high enough to overco me the interstellar plasma frequency 
l|Brodin et al.ll2001bf) . A numerical estimate for the case 
of a merging ns-ns binary shows that the amplitude of the 
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emws can be significant. Longitudina l waves are excited 
by higher-order GW-emw interactions l)Brodin et all2000|) 
and magnetos onic waves (msws) by GWs propagating in a 
low-/? plasma (jPaoadopoulos et"ai1l200l() . 

In this article we consider NS-NS mergers as a source for 
GRBs and apply the 3+1 formalism to the interaction of 
the Gws emitted by the merger with the ultra-relativistic 
wind of magnetized plasma around the binary. In the last 
seconds before the collapse to a black hole, a considerable 
fraction of M P) c 2 is released in to this plasma in the form 
of gws (see Ijanka et al.lEo02l Table 1). We show that 
these waves distort the extremely strong magnetic field, 
frozen into the plasma, and excite growing magnetosonic 
waves in the wind. Already before the merger, the binary is 
embedded in a relativistically expanding magnetized wind 
of (mainly) leptons from the orbiting neutronstars, so even 
a small transfer of GW energy to the wind may provide 
an interesting central engine mechanism to fuel a GRB 
fireball. 

The outline of this article is as follows. In Sect. the 
covariant expressions for the electromagnetic fields, the 
energy-momentum densities and the orthonormal tedrad 
for a linearly polarized GW are recapitulated. A closed set 
of linearized GRM equations in the metric of a GW is de- 
rived in Sect. Eland solved, first in the comoving frame of 
the plasma (Sect. 0} and then in the frame of an observer 
at rest with respect to the merger (Sect.[5J|. A numerical 
example and the interpretation of our results are given in 
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3. Magnetohydrodynamics in the comoving frame 



Sects. and [7| respectively. Sect. |H] comprises our conclu- 
sions. 

Throughout Secs.ErOO Gaussian geometrized units are 
adopted (c = 1) and Latin indices stand for 0,1,2,3. In 
Sec. El however, the numerical results are given in Si units. 

2. Covariant fluid equations 

2.1. Electromagnetic field equations 

Maxwell's equations in terms of the electromagnetic field 
tensors and the 4-current density = (r,j) are: 



ob 



47Tj° and V b T a = 0, (1) 

where the covariant Faraday tensor F ab and its dual !F ab , 
can be decomposed into 4- vectors that, in the rest frame of 
an observer with 4- velocity u a , reduce to the electric and 
magnetic field stren gths, E a = (0,E) and B a = (0, B) 
( Lic hnerowicall967j) : 
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In ideal mhd, the electric field vanishes in the rest frame 

in any 



of the plasma: E = 0. Therefore E a = F ab u b 
frame and Faraday's tensor reduces to: F ab 

2.2. Energy and momentum conservation 



-abed 



B c u d . 



Energy and momentum conservation follow from the 4- 
derivative of the energy-momentum tensor. For a magne- 
tized ideal plasma this can be expressed in terms of the 
proper matter energy density u, the pressure p and the 
metric g ab , as 



energy density /i, tnc pressure p and tJ 
llWeinherdll97ilHawking &r. Ellislll 973h : 



T ab = (fi + p)u a u b + pg ab + -t-(F a c F bc -\g ab F cd F cd ) 

An 4 

.b , 



V b T ab = V b [Qjl + p)u a u b + pg ab ] - F ab Jnlb = 0, (3) 
where the current density satisfies Eq. (TJJ. 



2.3. Tedrad system for a gravitational wave 

In the linearized theory of gravity one usually splits the 
metric: g ab = n ab + h ab , where r\ ab = diag(-l, 1, 1, 1) is the 
Minkowski metric, and h ab is the space-time perturbation 
caused by the GW. For a transverse-traceless, linearly (+) 
polarized, monochromatic GW with frequency w g = k g 
propagating in the z-direction, this field sa tisfies: h ab = 
du^ 0,/i,-/t,0), with h(z-t) = he^*'^ l(Misner et all 

The proper (observer) reference frame, however, co- 
moving with freely moving bodies is defined with respec t 
to the natural orthonormal tedrad l|Marklund et all20 0C0: 



e (0 f = (d u 0,0,0) , e (1 f 

e (2 f = (0, 0,(1 + |K, 0) , e (3 f 



(0,(1-1)^,0,0), 
(0,0,0,^). 



Decomposed with respect to this tedrad, the metric re- 
duces to that of flat space: g^ ab > = rf b and the GRM equa- 
tions closely resemble their Newtonian equivalents. 



The physical situation we want to consider is that of a 

perfectly conducting, ideal plasma in the presence of a 

background magnetic field along the x-axis, perpendicular 

to the direction of GW propagation (Fig.^l. First we study 

the plasma rest-frame where: B^ — Boe x , /i^ — p and 
w (o) =r fo) =p (o) =J -(o) = £(o) =a 

The effect of the GW is to induce small perturbations 
in all these quantities. Therefore, all equations will be lin- 
earized around the unperturbed state. 

3.1. Maxwell's equations 

The relevan t, linearized Maxwell e quations in the specified 
tedrad are l|Marklund et aljE00fj|) : 
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where the gravitationally induced current densities are just 
the collected Ricci rotation coefficients or GW-terms: 
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The electric field can be eliminated by assuming the ideal 
MHD approximation of a collisionless plasma (zero resis- 
tivity), where the electric field vanishes in the comoving 
frame and Ohm's law reduces to: 



-v^xB (0 \ 



(7) 



3.2. Conservation equations 



Charge continuity follows readily from the antisymmetry 
of F ab and Eq. Q: V a (V b F ah ) = iwVJ^ = 0. 

The evolution of the magnetic energy W and the 
Poynting flux S are just projections of Eqs. m^Jl onto 
B(°>: 



dz 



dt 

with F L =j x B the Lorcntz force, W (0) = B'fi 

W - ' 13 - 4i ' L — 3 Ta- 



rn 

(9) 
8ir and: 
x B<°>. 



Eq. (0 results in particle and momentum conservation in 
terms of the momentum, tv^ = /i^v^ and the spatial 
stress tensor T' 1 ': 



dt 



V 







9tt (1) 
dt 



V-T« = 0, 



(10) 



or, equivalently, as a linearized equation of motion (eom 
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dt 
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(11) 
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in terms of the sound velocity c 2 = T(p^ / p^), where 
r is the polytropic index. To first order, the EOM does 
not contain any GW terms and the coupling of the GW 
to the plasma occurs only through the Lorentz force as 
can be seen from Eqs. JH1 1 1 1 1) (or equivalently, through 
the current density which couples the EOM to Maxwell's 
equations I0H3J). 



4. Wave solutions in the comoving frame 

4.1. Laplace transforms 

The proper way to handle the unstable response of a 
plasma to a disturbance is to make use of Lapl ace trans- 
forms l|Landau fc Lifshit d Il975t lMelroselll986|) . By giv- 
ing the frequency or wavenumber a small positive imag- 
inary part, one allows for damping or growth and a way 
to deal with singularities due to resonance. In an initial 
value problem it is customary to Laplace transform with 
respect to time, whereas in a boundary value problem, as 
we are considering here, a Laplace transform with respect 
to space is more suitable. These are defined by: 



The inverse transformations lead to: 



*■(*) ^ So 



<f(z)dz , f(z) = ££l,e"F(8)d8, (12) 



where rj is an arbitrary positive constant chosen so that 
the contour of integration lies to the right side of all sin- 
gularities in F(s). Below, we will transform with respect 
to s = ik, where %i[k] w k g = u g , so the contour of 
integration follows the real axis except for the poles at 



fc = ±uj g /uA and k 



'9> 



where it deviates infinitesimally 



into the upper half imaginary plane. 

For the time dependence we use Fourier transforms, 
which implies that the perturbations oscillate at the fre- 
quency of the driving GW: oc exp (— iu> g t). 

4.2. Dispersion relations and wave solutions 

We will assume that the tenuous, strongly magnetized 
leptoid surrounding the binary has a low plasma-beta 
/3 P i = {cs/va) 2 = 47Tp/Sq <C 1, where the gas pressure 
is negligible with respect to the magnetic pressure. The 
classical Alfven speed v\ = B 2 /(Anp) > 1 so the dis- 
placement current is important. Perturbations propagate 
in the relativistic plasma with the generalized Alfven speed, 
= 1 + l/ui. 

In Laplace and Fourier space the 13 GRM equations 
derived in Sect. 13.21 reduce to an algebraic system with 5 
non-trivial solutions: 1 
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1 Since the only zeroth order quantities are So and p, the 
superscripts indicating the perturbations will be omitted. 
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Since Afc is very small (see Sect. we can expand the 
solutions around Afc = 0, or equivalently ma = 1, to find 
the dominant terms: 
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These are fast magneto- acoustic waves, propagating at the 
Alfven speed in the same direction as the GW and growing 
linearly with distance. 

5. Observer frame 

For a plasma flowing out relativistically in the z-direction 
with Lorentz factor 7 to t = — (f3+v z ) 2 ) ~ 7 + 7 3 /3v z 

(7 = jp) corresponding to a constant velocity (3 « 1, the 
full set of linearized GRM equations is: 



CONTINUITY 



7 3 p { m + Pik } v * = -iv B o 

2 §T-^T+ = ^(B +E )e^-^ AMPERE (17) 



dE v 
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E 
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where all quantities are now defined in the observer frame 
at rest with respect to the binary (primes are used to 
indicate comoving quantities only where there is the risk 
of confusion). Explicitly: p, (3, Bo and Eq = E y °^ are the 
zeroth order quantities and p, v z , B x , E y , j y , h are the 
perturbations. The most important difference with respect 
to the comoving frame is the factor (1 — (3)Bq in the GW 
terms, due to the background electric field seen by the 
observer. 

The solutions to the equilibrium system without the 
GW-terms, are again oc e luJ g( z / u A-t) UA _ the 
Lorentz boosted, generalized Alfven speed. The full sys- 
tem can be solved along the same lines as in Sect. 14.21 
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by straightforward calculation using Laplace and Fourier 
transformations. The physical differences between the co- 
moving and the observer frames are clear from Eqs. I|17|l . 

A different approach - leading to the same result - is 
to Lorentz transform Eqs. (|14|) . Since Lorentz transforma- 
tions are linear transformations and furthermore the phase 
of a plane wave is an invariant, the general solutions have 
the same form as Eqs. Ijf4|l . but with different amplitudes. 

The msw components excited by a GW propagating 
through a rclativistically flowing plasma are given in terms 
of observer quantities by: 
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where the constants A-Ag,A are defined in Appendix 1X1 
As in Sect. 14.21 the dominant terms are found by expand- 
ing around Ak = (justified in Sect. 16. 3|) : 
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where: p = " A . Apparently, the interaction is less 

efficient when the plasma is escaping relativistically, than 
when it is at rest with respect to the source of Gws. 

As a final remark: Eq. I|19[) is equivalent to the result 
of an initial value approach, where the Laplace transfor- 
mations are performed with respect to t instead of z. The 
amplitudes in Eq. 1|19[) become proportional to tot instead 
of kz, but the characteristic timescale T is related to the 
size of the interaction region R by: loT ~ u(R/ua) = kR. 

6. Numerical estimates 

In this section the results obtained in the previous sections 
are applied to the environs of a NS-NS binary close to merg- 
ing. Both stars will have an electron-positron wind filling 
the surrounding space with plasma up to large distances 
(whether the plasma is e ± or baryon loaded is not impor- 
tant as long as it satisfies the ideal mhd condition). The 
extent of the interaction region is determined by the dis- 
tance i? max from the source where either the force-free or 
the z«l/ Afc assumption brakes down. Within this scale 
height, we estimate the magnitude of the excited msw am- 
plitudes, including the decrease of the GW-amplitude, the 
density and the magnetic field with distance. A short nu- 
merical analysis is made for the magnetic fields, plasma 
densities, Lorentz factors and Alfven velocities in the rel- 
ativistic plasma wind, summarized in Tabled 



Note: In this section proper dimensions in c are restored 
and the numerical results are converted to Si units. 



6.1. Magnetic field configuration 

For pulsar-like neutronstars, the magnetic field close to 
the surface falls of as a dipole: B(r) = B+ (R*/r) 3 . In the 
mhd approximation the plasma is 'glued' to the field lines 
and forced to corotate up to the lightcylinder where coro- 
tation requires superluminal velocities. Here, we consider 
a NS-NS merger where each NS has its own magnetic field. 
The key element in the electromagnetic description of ro- 
tating magnetized stars is the deviation from inertial mo- 
tion. In our case this is not the usual stellar rotation but 
the orbital motion combined with the individual stellar ro- 
tations. As the merger coalesces the orbital frequency in- 
creases and dominates over any other (rotational) motion. 
Therefore, we assume that at the end of the spiral- in phase 
the orbital rotation of the binary (with flh ~ I0 3 rad/s) de- 
termines the light-cylinder radius: R\ c = c/fib — 300km. 
Here the field lines, anchored on the stellar polar caps, 
open up and the plasma is free to flow out along the field 
in a force-free wind in which the toroidal component of th e 
field dominates the poloidal component l|Kuiipersl [200f ): 



(20) 



B t (r>R ic ) - B lc (^) , B p (r>R lc ) = B lc (^) 2 , 
Sic = B,(^) 3 ~ 4 • I0 5 T(^, 9 )(O b , 3 ) 3 (i?,,4) 3 . 

6.2. Nature of the wind 



Above the polar caps of a pulsar the field lines are open, 
particles flow outwards and a steady charge density can- 
not be maintained at the Goldreich- Julian density every- 
wher e no j = f2b--B/(27rec) ijBeskin et al Jl993HLvubarskvl 
1995). As a result, a strong electric field develops along the 
magnetic field above the polar cap and charged 'primary' 
particles are extracted from the surface with a density 
n p ~ ^gj* and accelerated to high Lorentz factors (a typ- 
ical number is 7 p ~ 1 7 ). Note that the available potential 
jump is proportional to A4" oc fib-B* in our case and can 
be much larger than for single pulsars. 

Processes such as curvature radiation and inverse 
Compton emission then result in a cascade of 'secondary' 
e pairs with a particle number density n s = Mn p , where 
M is called the multiplicity. Due to energy conservation 
n p 7p = n s 7 s , so the Lorentz factor of the secondary par- 
ticles is 7 S = 7 P /M ~ 100 for M ~ I0 5 . The secondary 
plasma flows out as a relativistic wind along the open 
magnetic field lines which have a dipole-likc radial de- 
pendence up to R\c, and develops into a spiral (Eq. I20f> 
further out. The charge density in the wind is adjusted 
to the local Goldreich- Julian density everywhere and the 
wind remains force-free up to a large distance. 

At the lightcylinder the density is: 



Sic 



2-Kec 



4- I0 2 



(21) 
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Lab/Observer frame 
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Plasma/Comoving frame 



fe g c 
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2c n. 



fcg ~ Hr 7 m -1 
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Table 1. Formulas and numerical values in the frames at rest with respect to an observer and the plasma, respectively. 
In the numerical examples of this table we have taken: £?* = 10 n T, 7 = 7 S = 100, M = 10 5 , n^/2TT = 160Hz. 



6.3. Alfven speeds &. Interaction lengthscale 



we find: 



The classical Alfven speed va is proportional to Bq 
In the force free wind we have Bq oc l/r and n cx l/r 2 , 
so va(t > R\c) — constant. Therefore, we can evaluate 
the Alfven velocities at the light-cylinder in the comoving 
frame (using Eqs. JUJ): 



Bi 



^ and 
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Aim', m c c 2 



7s 



B. 



a, 



47T7 s ni c m e c 2 



27 p f2 b 



»1, 



(22) 



in terms of the electron cyclotron frequency O e = . For 
the generalized Alfven speed and the wavenumber differ- 
ence Ak' , one finds (see also Table QJ: 



Ak' = k' - ki = 



1 



< 1. 



»1, 



(23) 



The equivalent quantities in the observer frame are de- 
rived from Lorentz transformations and can be found in 
Table ^ The important conclusion from these estimates 
is that the maximum distance over which the MSW am- 
plitudes can grow is: z <C 1/Afc ~ 10 15 m (see Eq. H19[l). 
which is comparable to the limits on the MHD approxima- 
tion as calculated bv lSpruit et"al] l|200ljL 

6.4. Magnitude of excited magnetosonic waves 

The frequency of Gws emitted by a merging binary is re- 
lated to the angular f requency of the last orbits by: t-u g = 
2f7 b -47r-10 3 rad/s l)Shibata fc Urvul2002|) . so ii lc ~50km 
and i?i c ~10 7 T. We assume that for i?; n ~10i?i c we are in 
the far field of the merger and estimate the amplitude of 
B^ (z, t) in Eq. I|19(l at the start of the interaction region. 
Taking into account the (approximate) l/r decrease of h 
and Bq, 



h ~ 10~ 3 (^) , B ~ Bu 



B^(r>R in ) ~ 0.5T(Sa) 



r u 



(24) 



The volume integrated energy of the MSW grows linearly 
with distance since B { V ocB oc l/r and V cx r . For an 
interaction region of -R max ~ 0.03 pc we find for the mag- 
netic component of the excited MHD wave a total energy 
of: 



T (l) _ V B„B^ 



10 dY J, 



(25) 



which amounts to a fraction 10 6 of the magnetic energy 
Tg = ~ 10 43 J of the Poynting flux dominated wind 
in the same volume. 




Figure 1. A GW propagating in the positive ^-direction 
across an ambient magnetic field (in the x-direction) ex- 
cites a MSW. The orientations of the MSW components are 
indicated in the inset. 
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7. Discussion 

It is known that, to first order in h, GWs propagating along 
a uniform magn etic field do not couple to the field, nei- 
ther in vacuum ( Boccaletti et al . 1970), nor in a plasma 
l|Brodin et al.l2000HPapadopoulos et al.l200l[) so we stud- 
ied GW-propagation perpendicular to the magnetic field 
frozen into a plasma. An intuitive illustration of the pro- 
cess is given in Fig.^ where a GW traveling in the positive 
z-direction, deforms (imaginary) rings of plasma-particles 
in the x-y plane into ellipses with long axes alternating 
periodically along the magnetic field direction (x-axis) and 
the y-axis. Consequently, the uniform field is periodically 
compressed and stretched leading to a modulated magnetic 
field strength even though the GW is divergence-free. The 
resulting magnetic pressure gradients in the z-direction 
try to re-establish a uniform field configuration and ex- 
cite the compressional matter part of the MSW. Since the 
plasma is glued to the field lines it is dragged along with 
the field. The velocity in the ^-direction generates an elec- 
tric field (—v x B) with a corresponding E x B drift- 
velocity, whereas the magnetic gradient induces a B x VI? 
current density and Lorentz force. 

The problem o f GW propagation in a pla sma at rest 
was also studied bv lPapadopoulos et al.l l)200l|) who found 
that the excitation of msws by a GW is only possible when 
there is a wavenumber mismatch: B% oc Afc. This does 
not agree wit h our results (as presen ted in Sec. 14.2(1 . nor 
with those of iMarklund et alJ l)2000[) who show that the 
growth rate i n a tenuous plasma sm oothly matches that 
in a vacuum l(Boccaletti et al.|ll97fllk Both from the GRM 
equations (Eqs. or fTTjl ) and from their solutions 

(Eqs. Ijl 60 or G29) it is clear that in the vacuum limit 



fi {1) tx jy oc p I vanish as does v\ L> . The GWs excite 
EMWs that propagate with the speed of light (Afc J. 0) and 
with amplitudes that can be obtained from Eq. iffifll or 
Eq. (|19jl by taking the limit ua — > c and k — » fc g (and 
7 = 1). This result is not surprising as we have taken 
into account both the material current and the vacuum 
displacement current. 



(i) 



8. Conclusion 

The results presented in this paper are a first attempt to 
study the interaction of the relatively strong GWs emitted 
by a GRB with an ultra-relativistic plasma wind. The space 
surrounding a merging NS-NS binary is already filled with 
such a wind up to large distances (~ 0.1 pc). Moreover, 
in the merger almost all of the binding energy is released 
in the form of GWs. 

We derive a closed set of GRM equations both in the 
natural orthonormal measurement frame (the 3+1 split) of 
the plasma and for an observer at rest with respect to the 
binary. These non-coordinate equations strongly resemble 
their Newtonian equivalents but have extra source terms 
due to the GW. These gravity terms act as a driver for fast 
magnetosonic waves with amplitudes that grow linearly 
with distance and are proportional to the GW-frequency 



and amplitude and to the ambient magnetic field strength. 
It is the extended force-free wind in which the magnetic 
field only falls of as 1 jr that provides the long interaction 
lcngthscale. 

The total amount of energy that is transferred from the 
GWs to the plasma, as given in Eq. I|25[) , is substantial but 
for this case still much smaller than the average observed 
GRB energy. Note that for magnetars one can have ~ 
10 12 T and lu s ~ 15 kHz, so that [BqB^/Atx could be as 
much as a factor 10 7 larger. However, for magnetars it is 
not obvious what to assume for the surrounding plasma. 

In future work we will investigate what the observable 
effects of the GW-MSW interaction are on the emitted radi- 
ation and its polarization and how a GW- chirp alters the 
results. 
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Appendix A: Constants in Eq. ( )18( > 

a h l+u A (1-/3) 2 ^ h y g 

ly — 4 1— u A 1 — u A /3 — 4 7 2 Afe 

Ai = 2 ^^ (i ' 4 " A :r ( \tiv (1 " ffl) - i+^[Afc] 2 

a _ q 1-u a (1+« a +2/3 2 -(1 + (4- ma ) ma )/3 _ 1 /-or A il 2 
ii3 _ Z (l-/3) 3 (l+u A ) 2 — L-\-U[lSK\ 

i T ^|±f) 2 (/3 2 -2u A /?+l)=0[Afc] 2 

a 5 = 4 (1 ^/ )a ;r-^ = i-o[^] a 

Ae = (^^f) 2 =0[Afc] 2 

a 7 = p&p l x u _f = i + om 

a 8 = (m) 2 ^& =o[A k r 

A 9 = 4 7 2 (^_£) 2 = l + G[Afc] 
A = ^T+J =7 2 (1-^A/3)A 



A4 ~ ^ V !+«A 1-/3 J 



H 



ua ua(1+P )— 2/3 u a 1 
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